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Abstract 

We prove that the diffusion coefficient for the two dimensional asymmetric simple 
£Sj \ exclusion process diverges as (log£) 2 / 3 to the leading order. The method applies to 

nearest and non-nearest neighbor asymmetric simple exclusion processes. 

in 

1 Introduction 

o 

Asymmetric simple exclusion process is a Markov process on {0, 1} Z " with asymmetric jump 
rates. There is at most one particle allowed per site and thus the name exclusion. The 
Q-r particle at a site x waits for an exponential time and then jump to y with rate p(x — y) 
provided that the site is not occupied. Otherwise the jump is suppressed and the process 
starts again. The jump rate is assumed to be asymmetric so that in general there is net 
drift of the system. The simplicity of the model has made it the default stochastic model 
for transport phenomena. Furthermore, it is also a basic component for models [Q [0 with 
incompressible Navier-Stokes equations as the hydrodynamical equation. 

The hydrodynamical limit of the asymmetric simple exclusion process was proved by 
Rezakhanlou Jl3| to be a viscousless Burgers equation in the Euler scaling limit. If the 
system is in equilibrium, the Burgers equation is trivial and the system moves with a uniform 
velicity. This unifrom velocity can be removed and the viscosity of the system, or the 
diffusion coefficient, can be defined via the standard mean square displacement. Although the 
diffusion coefficient is expected to be finite for dimension d > 2, a rigorous proof was obtained 
only a few years ago by estimating the corresponding resolvent equation. Based on the 
mode coupling theory, Beijeren, Kutner and Spohn conjectured that D{t) ~ (logt) 2 / 3 
in dimension d = 2 and D(t) ~ t 1//3 in d = 1. The conjecture at d = 1 was also made by 
Kardar-Parisi- Zhang via the KPZ equation. 

This problem has received much attention recently in the context of integrable systems. 
The main quantity analyzed is fluctuation of the current across the origin in d — 1 with the 
jump restricted to the nearest right site, the totally asymmetric simple exclusion process 
(TASEP). Consider the special configuration that all sites to the left of the origin were 
occupied while all sites to the right of the origin were empty. Johansson H observed that 
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the current across the origin with this special initial data can be mapped into a last passage 
percolation problem. By analyzing resulting percolation problem asymptotically in the limit 
iV — > oo, Johansson proved that the variance of the current is of order t 2 / 3 . In the case of 
discrete time, Baik and Rains analyze an extended version of the last passage percolation 
problem and obtain fluctuations of order t a , where a = 1/3 or a = 1/2 depending on 
the parameters of the model. Both the approaches of || and are related to the earlier 
results of Baik-Deift- Johansson [0] on the distribution of the length of the longest increasing 
subsequence in random permutations. 

In |10j (see also ||11||), Prahofer and Spohn succeeded in mapping the current of the 



TASEP into a last passage percolation problem for a general class of initial data, including 
the equilibrium case considered in this article. For the discrete time case, the extended 
problem is closely related to the work 0, but the boundary conditions are different. For 
continuous time, besides the boundary condition issue, one has to extend the result of 
from the geometric to the exponential distribution. 

To relate these results to our problem, we consider the asymmetric simple exclusion 
process in equilibrium with a Bernoulli product measure of density p as the invariant measure. 
Define the time dependent correlation function in equilibrium by 

S(x,t) = 7^,(0)) 

We shall choose p = 1/2 so that there is not net global drift, ^2 X xS(x, t) = 0. Otherwise a 
subtraction of the drift should be performed. The diffusion coefficient we consider is (up to 
a constant) the second moment of S(x,t): 



X 

for large t. On the other hand the variance of the current across the origin is proportional 
to 

^2\x\S(x,t) (1.1) 

x 

Therefore, Johansson's result on the variance of the current can be interpreted as the spread- 
ing of S(x, t) being of order t 2 ^ 3 . The result of Johansson is for special initial data and does 
not directly apply to the equilibrium case. If we combine the work of |1(| and 0, neglect 
various issues discussed above, and extrapolate to the second moment, we obtain a growth 
of the second moment as t 4//3 , consistent with the conjectured D{t) ~ t 1 / 3 . 

We remark that the results based on integrable systems are not just for the variance of 
the current across the origin, but also for its full limiting distribution. The main restrictions 
appear to be the rigid requirements of the fine details of the dynamics and the initial data. 
Furthermore, it is not clear whether the analysis on the current across the origin can be 
extended to the diffusivity. In particular, the divergence of D{t) as t — > oo in d = 1 has not 
been proved via this approach even for the TASEP. 

Recent work of || has taken a completely different approach. It is based on the analysis 
of the Green function of the dynamics. One first uses the duality to map the resolvent 
equation into a system of infinitely-coupled equations. The hard core condition was then 
removed by using variational principles. Once the hard core condition was removed, the 
Fourier transform becomes a very useful tool and the Green function was estimated to degree 
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three. This yields a lower bound to the full Green function via a monotonicity inequality. 
Thus one obtains the lower bounds D(t) > i 1 / 4 in d = 1 and D(t) > (logt) 1 / 2 in d = 2 [|J. 
In this article, we shall estimate the Green function to degrees high enough to determine the 
leading order behavior D(t) ~ (log t) 2 / 3 in d = 2. 



1.1 Definitions of the Models 

Denote the configuration by rj = {i] x ) xe x d where rj x — 1 if the site x is occupied and r] x = 
otherwise. Denote r/ 1 ^ the configuration obtained from rj by exchanging the occupation 
variables at x and y : 

{r/ 2 if z ^ x, y, 
^ if z = y and 
77^ if £ = x. 

Then the generator of the asymmetric simple exclusion process is given by 

d 

(Cf)( v ) = ^^K^z/Kti-^t/^)-/^)]- (i.2) 

3=1 x&Z d 

where {ek, 1 < k < d} stands for the canonical basis of Z d . For each p in [0, 1], denote by u p 
the Bernoulli product measure on {0, l} zd with density p and by < •, • > p the inner product 
in L 2 (i> p ). The probability measures v p are invariant for the asymmetric simple exclusion 
process. 

For two cylinder functions /, g and a density p, denote by (/; g) p the covariance of / and 
g with respect to v p : 

(f;g) P = (fg)p- (f) P (g) P - 

Let P p denote the law of the asymmetric simple exclusion process starting from the equilib- 
rium measure v p . Expectations with respect to P p is denoted by E p . Let 

S p (x,t) = E p [{ Vx (t) - %(0)}%(0)] = (vMVo(0))p 

denote for the time dependent correlation functions in equilibrium with density p. The 
compressibility 

X = X(p) = ^2(Vx] Vo) P = Yl *) 

X X 

is time independent and xip) = P(l — P) lY1 our setting. 

The bulk diffusion coefficient is the variance of the position w.r.t. the probability measure 
S p (x, t)x _1 m ^ d divided by t, i.e., 

D ij(p,t) = t| XiXjS^x^x" 1 - (vit)(vjt)\, (1.3) 
^ xez d ' 

where v in M. d is the velocity defined by 

vt = YxS^x,^ 1 ■ (1.4) 

x&L d 
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For simplicity, we shall restrict ourselves to the case that the jump is symmetric in the 
y axis but totally asymmetric in the x— axis, i.e., only jump to the right is allowed in the 
x— axis. Our results hold for other jump rates as well. The generator of this process is given 
by 
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(Cf)( v ) = ^Iv^-v.+eJifiv^-fivV + ^ifiv^-fiv)) 



(1.5) 



where we have combined the symmetric jump in the ?/-axis into the last term. We emphasize 
that the result and method in this paper apply to all asymmetric simple exclusion processes; 
the special choice is made to simplify the notation. The velocity of the totally asymmetric 
simple exclusion process is explicitly computed as v — 2(1 — 2p)ei. We further assume that 
the density is 1/2 so that the velocity is zero for simplicity. 

Denote the instantaneous currents (i.e., the difference between the rate at which a particle 
jumps from i to i + e ; and the rate at which a particle jumps from x + ti to x) by w x>x+ei '- 



Tj x [l - Tj x+ei ] 



Vx+e2 Vx 



1.6) 



We have the conservation law 



i=l 



Let Wi{rf) denote the renormalized current in the i-th direction: 

Wiiv) = w , ei ~ (wo, ei ) P ~ jE(w , ei )e {vo ~ p) 

do o=p 

Notice the subtraction of the linear term in this definition. We have 

wi{ri) = (r) - p){r) ei - p) + p[v ei - Vo], w 2 (r)) = 

Define the semi-inner product 

((g, h )) P = < T x9 ; h > p = Y < T * h ; 9 > P ■ 

x& d x&Z d 



;i.7) 



Since the subscript p is fixed to be 1/2 in this paper, we shall drop it. All but a finite number 
of terms in this sum vanish because v p is a product measure and g, h are mean zero. From 
this inner product, we define the norm: 



'U 



Notice that all degree one functions vanish in this norm and we shall identify the currents 
w with their degree two parts. Therefore, for the rest of this paper, we shall put 



= (Vo - p){Vei - p), w 2 (rj) = 



1.9) 



Fix a unit vector £ e Z d . From some simple calculation using the Ito's formula we 
can rewrite the diffusivity as 

2 



2 X 



r 1 ' 2 / da • w){r,(a)) 



1.10) 
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r 1/2 / dswMs)) 



This is some variant of the Green-Kubo formula. Since W2 = 0, D is a matrix with all entries 
zero except 

1 1 

Du = - + - 

2 X 

Recall that f Q °° e~ xt f(t)dt ~ X~ a as A — > means (in some weak sense) that /(£) ~ 
Throughout the following A will always be a positive real number. The main result of this 
article is the following Theorem. 

Theorem 1.1 Suppose that the density p = 1/2. Then there exists a constant 7 > so that 
for sufficiently small A > 0, 

POO 

\- 2 \ log A| 2 /V^ losloslosA l 2 < / e- xt W u {t)dt < \- 2 \ log A| 2 /Vl losloglogA l 2 

Jo 



From the definition, we can rewrite the diffusion coefficient as 

W n (t) = I + - I [ ((e u V, w 1 )) duds 
1 X Jo Jo 

Thus 



co pt ps 

e~ M W n (t)dt = -— + -! dt I I e~ xt ({e^w^wi}) duds 



2A 2 x Jo Jo Jo 

I ^ /* oo s poo 



n H — / du 
2A 2 x Jo 

1 

2~A 



dt e - x(t - u) (J ds ) } ((e- Au e"V , ^x)) 

(1.11) 



Therefore, Theorem |LT] follows from the following estimate on the resolvent. 
Theorem 1.2 There exist a constant 7 > such that for sufficiently small A > 0, 
|logA| 2/3 e- 7|losloglosA|2 < ((wt, (A - Cy 1 w 1 }} < |i gA| 2 /V |losloslogA|2 



From the following well-known Lemma, the upper bound holds without the time integra- 
tion. For a proof, see ||. 



Lemma 1.1 Suppose /i is an invariant measure of a process with generator C. Then 



r 1 ' 2 w(rj(s))ds) J < ((^(r 1 -£)->) 



;i.i2) 



Since ui\ is the only non- vanishing current, we shall drop the subscript 1. 
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2 Duality and Removal of the Hard Core Condition 



Denote by C = C(p) the space of z/ p -mean zero cylinder functions. For a finite subset A of 
Z d , denote by £a the mean zero cylinder function defined by 

Denote by Ai n the space of cylinder homogeneous functions of degree n, i.e., the space 
generated by all homogeneous monomials of degree n : 

M n = \h e C ; h = , h A e r} . 

|A|=n 

Notice that in this definition all but a finite number of coefficients h\ vanish because h is 
assumed to be a cylinder function. Denote by C n = Ui<j< n A4j the space of cylinder functions 
of degree less than or equal to n. All mean zero cylinder functions h can be decomposed 
as a finite linear combination of cylinder functions of finite degree : C = U n >iA4 n . Let 
£ = S + A where S is the symmetric part and A is the asymmetric part. Fix a function g in 
■M-n : g = J2 A |A|= n ^A^A- A simple computation shows that the symmetric part is given by 

(Sg)(v) = ~2 X/ X {gnu{x+ ei } - gnu{x}} &2u{x+ei} - £nu{x} 
x&L d n,\n\=n-i 

Qn{x,x+ei}=<p 

The asymmetric part A is decomposed into two pieces A = M + J so that M maps Ai n into 
itself and J = J + + J_ maps M. n into A4 n ~i U M. n +i- 

(m 9 )w = ^ E £ { 



xez d n,|n|=n-i 

Hn{a;,a;+ei}=0 



(J+gXv) = ~VP( 1 - P) ^ ^ |^U{x+ei} - ^nu{x}}^nu{x,x+ei} , 



xez d n, |n|=n-i 

f!n{x,x+ei}=<j!> 

(J-g)(v) = -Vpi 1 - P)Y1 Yl {gnu{x+ ei } -gnu{x}}tn ■ 

xe1 d Q, |Q|=n-l 
f2n{x,x+ei}=0 

Restricting to the case p = 1/2, we have M = and thus J = A. Furthermore, 
A* + = —A_. We shall now identify monomials of degree n with symmetric functions of n 
variables. Let E\ denote the set with no double sites, i.e., 

£i = {x„ := (xi, • • • , x n ) : Xi ^ x 5 , for i ^ j} 

Define 

f(x u ■■■ ,x n ) = f{ xlt ..., Xn y, if x n G 

= 0, XXn?£l. (2.1) 
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Notice that 



'{[E/^a] 3 } = ^ E i/(*i,-,*»)i s 



^ , , 

n! 

|A|=n xi, ■■■ ,x n £^ 



From now on, we shall refer to f(x\, • • • , x n ) as homogeneous function of degree n vanishing 
on E\. 

With this identification, we have 

u'i(0,ei) = i/;i(ei,0) = 1 

and zero otherwise. Since we only have one non- vanishing current, we shall drop the subscript 
1 for the rest of this paper. 

If g is a symmetric homogeneous function of degree n, we can check that 



A + g(x u ■ ■ ■ ,x n+1 ) 

^ n+1 

— — ^ ^ ^ b(^l) ■ ■ ■ j x i e l> ■ ■ ■ j ^n+l) 9{ X 1j > ^i) ' ' ' ) 2" J) 5 2-n+l)] 



i=l j^i 

x^Xj - a?i - ei) JJ (1 - 5(a?j - x fc ) ) (2.2) 
where 5(0) = 1 and zero otherwise. We can check that 

n 

Sg(xt, ■■■ , x n ) = a 2J ]^[ ( 1 - + °" e a - ^a) ) 

i=l ct=± a=l,2 fc^j 

x • • ■ Xi + ae a , •■•,£„)- • • ■ , £„)] (2.3) 

where a is some constant and 5(0) = 1 and zero otherwise. The constant a is not important 
in this paper and we shall fix it so that S is the same as the discrete Laplacian with Neumann 
boundary condition on S\. 

The hard core condition makes various computation very complicated. In particular, the 
Fourier transform is difficult to apply. However, if we are interested only in the orders of 
magnitude, this condition was removed in ||. We now summarize the main result in ||. 

For a function F, we shall use the same symbol (F) to denote the expectation 



i E F(xi,---,x n ) 



n. 

XI,--- ,x n £l? 

We now define A + F using the same formula except we drop the last delta function, i.e, 

n+1 

A + F(x ir -- ,x n+ i) = --^2^2^F(x ir -- ,Xi + e ir -- ,x jr --x n+1 ) 

i=l jyt% 

-F(x u ■ • ■ ,x h ■ • ■ ,xj, • ■ ■ ,x n+1 ) 5(x j -x i -ei) (2.4) 



Notice that (A + F) = 0. Thus the counting measure is invariant and we define A_ = — A* + , 
i.e., 

(A_G,F) = -(G,A+F) (2.5) 
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Finally, we define 

L = A + A, A = A + + A_, 
where the discrete Laplacian is given by 

n 

. Xi -\- <7€ a , . . . , X n ) F{Xi, . . . , Xj, , . . . , X n )] 

i=l cr=± o=l,2 

For the rest of this paper, we shall only work with F and L. So all functions are defined 
everywhere and L has no hard core condition. 

Denote by 7r„ the projection onto functions with degrees less than or equal to n. Let L n 
be the projection of L onto the image of ir n , i.e., L = 7T n Ln n . The key result of || is the 
following Lemma. 



Lemma 2.1 For any A > fixed, we have for k > 1, 

C-'k-^dw^L-^w}) < ((w^-'w)) < Ck 5 ((w 7 L-£w) 



(2.6) 



The expression ((w, L^w)) was also calculated in The resolvent equation (A- 
w can be written as 

(A - S)u n - A + tt n _i =0, 
A+Uk+i + (A — S)uk — A + Uk-i =0, n — 1 > k > 3 
A* + u 3 + (A - S)u 2 =w. 

We can solve the first equation of (|2.7|) by 

u n = (A - S')~ 1 A h ii n _i. 

Substituting this into the equation of degree n — 1, we have 



■L n )u 



(2.7) 



Solving iteratively we arrive at 

U2 



(X-S) + AU\-S)- l A, 



-1 



U n -2- 



(\-S) + A* + {(\-S) + --- + A* + ((\-S) + A* + (\-S)- 1 A + ^ A + } A, 
This gives an explicit expression for ((w, (A — L 3 ) _1 U7)), for example, 
((w, (A - L 3 )~ l w)) = ((w, \X-S + A* + (X- Sy'A, 



-1 



w . 



-1 



w >>. 



\-S + AU\-S + A*J\-S)- 1 A+\ A 



-1 



-1 



w 



(w, (A - L A )- l w)) 
(w, (A - L,y l w)) 

= ((w, \\ - S + A* + \ A - S + A; [A - S + A;(A - S , )- 1 A+]- 1 A t . !> "w)). 



(2- 



-1 



-1 



This is the expression we shall estimate for the rest of the paper. 
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3 Main Estimate 



We now introduce the following convention: Denote the component of p by (r, s). Denote 
Pn = {pir ■ ■ ,Pn),r n = {n, - ■ ■ , r n ) and s n = (si, • • • , s n ). The Fourier transform of 

[F(xi + ei, • • • , x n ) - F(x u • • • , x n )]S(x n+1 = x 1 + e x ) 

is given by 

^[F{ Xl + ei, • • • , X n ) - F(Xl, ■■■ , Xn )] e -^Pl+-+*nPn + (xi+ei)Pn + l] 
x 

= [e iri - e~ irn+1 ] F(pi + p n+1 , ■ ■ ■ , p n ) ~ [i(n + r n+ i)] F( Pl + p n+1 , • • • , p n ) 

All functions considered for the rest of this paper are symmetric periodic functions of period 
2tt. 

Since F is symmetric, we have 

71+ 1 

A^(Pn+l) = -J^^ -e~ lrm )F{ Pl ,--- ,Pj+Pm,--- ,PZ,--- ,Pn+l) (3.1) 

j<m 

We can also compute the discrete Laplacian acting on F 

n 

^F(Pn) = - E E ~ 2 + P ^ * * * = ^(P-)^(Pn) 

3=1 k=l,2 

where u?(p n ) = YJj=i w {Pj) and 

wfo) = - [ e ir * - 2 + e-^] - [e^ - 2 + e"*'] (3.2) 
We shall abuse the notation a bit by denoting also 

u(rj) = - [e ir * - 2 + e-^] 
Notice that yju{x) — iff x = mod 7r. When x ~ mod 7r, we have 

u{x) ~ | sina;| (3.3) 

By definition, we have 

{(F,G)) = ^2— E p ( x ir-- ,x n )G( Xl + z,--- ,x n + z) 



ni 

3*1 )*** l^TJ 



= E^ / d Pl ---dp n F( Pl ,--. ,p n )G( Pl ,--- ,p n )e^ 1+ - + ^ 

= ~\ J dp 1 ---dp n 5(p 1 ^ hp„)F(pi, • • • ,p„)G(pi, • • • ,p n ) (3.4) 

In other words, when consider the inner product ((•,•)), we can consider the class of 
F{pi,--- ,p n ) defined only on the subspace J2jPj = m °d 27r. We shall simply use the 
notation ^ . pj = to denote the last condition. 

From now on, we work only on the moment space and all functions are defined in terms 
of the momentum variables. Let d/i n +i(Pn+i) denote the measure 

n+l Ti+l 

dA*n+l(Pn+l) = 7- x x, ^ (E^) II ( 3 - 5 ) 
^ J"=l 3=1 
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3.1 Statement of the Main Estimate 

Let r be a positive constant and define 

^(p„) = {u,(p„)<|logA|- 2 -} (3.6) 
Denote the complement of Q T by B T . Define for k > the two operators 

Vl T ( Pn ) = «(r B )| log(A + w(p n ))r, Pn G £ r ; 

= w(r n )|log(A + w(p n ))|, Pn eB T (3.7) 

V^ T ( Pn ) = ui(r B )| log(A + u>(pn))\ K , Pn e G T ; 

= -|loglogA| 2 u,(r n ), p n G B T . (3.8) 

The main estimates of this paper are the following Theorem. 
Theorem 3.1 Let k and r be nonnegative numbers satisfying 

< k < 1 < t (3.9) 
Let n be any positive integers such that 

n 10 < | log log A | 1/2 (3.10) 



Suppose that for some 7 < | log log A | 



-3 



^« + i>7VS (3.11) 

as an operator. Let 

k = 1 - k/2 . (3.12) 

Then 

A* + (\ - S n+1 + n n+1 y 1 A + < 7 - x | loglog A| 2 U? T (3.13) 

as an operator. 

On the other hand, if 

< 7 _1 US' ( 3 - 14 ) 

Then we have 

A* + (X - S n+1 + > C 7 V? 2r (3.15) 

as an operator. 
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4 Upper Bound 

We first recall that for any two positive operators A, B, we have 

< A < B iff < B^ 1 < A^ 1 . 
Furthermore, the map B — > C*BC is monotonic. For 7 < | log log A|~ 3 , we have 

w(p B+1 )+ 7 VS(p n+ i)>0 



Thus we can substitute the Vt in Theorem p\l| by either V or U in the proof. 
By definition, 

{{F,A\\-S n+1+lY ^ r x AF))= f d^ n+1 ( Pn+1 ) JY^l'^f , 

J A + w(p n+ i) + 7V K + T ( Pn+1 ) 

Let V^' 2 ™ +1 denote the positive and negative parts of V™ Then 

A + W (p n+1 )+7VS(p n+1 ) > (l- 7 )^(p„ + i)+7V:;™ +1 (p n+1 )+ 7 [A + a;(p n+1 ) + V^(p n+1 )] 
Since 7 < | log log A |~ 3 , we have 

(1 - 7 )«(lWi) + 7V-'™ +1 (Pn + i) > . 

Thus we have 

((F, A*(A-5 n+1 + 7V^)- 1 AF)) 

< 7 / cf// n+1 (p n+1 ) — — . (4.1) 

J A + w(p n+1 ) + V K ; 2r ( Pn+ i) 

We now divide the integration into the good region Q 2r (p n+ i) and the bad region i3 2r (p rt+ i) 
In the good region, 

V+£ +1 (p n+ i) =uKr n+1 )|log(A + u;(p„ +1 ))r • 
Thus the contribution is 

[ A , \ /">2t / \ \A+F(p U --- ,p n +l)\ 2 (A „s 

/ d/x n+ i(p n+ i) £ Pn+i x , / \ 1 / vn 7T , — 7 7yn: • 4 - 2 

J A + u>(p n+1 ) + w(r n+ i)| log(A + w(p rt+ i))| K 

Since V^^ 1 = in the bad region, the contribution from this region is 

a 1 \Y2iri \ \A+F(pi, ■ ■ ■ ,p n+ i)\ 2 

d/x„+i(p„+i) B (p n+ i) — r 4.3 

A + oj(p n+ i) 
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4.1 Decomposition into Diagonal and Off-Diagonal Terms 

Denote by Q K the function 

©«(Pn+i) = [A + w(p n+ i) + w(r n+ i)| log(A + u:(p n+1 ))\ K }- 1 (4.4) 

The contribution from the good region can be decomposed into diagonal and off-diagonal 
terms: 



J d//„+i(p„+i) g 2r (p n+1 )Q K (p n+1 )\A + F(p 1 , ■ ■ ■ ,p n+1 ) 



(4.5) 



+ n(n-l)(n + l)«F, ^ F)) + ^ = ^ 2 * n + (( F, *^F» 
where 

«F, /C^F)) = | ^ n+1 ( Pn+1 ) K ( Pn+1 )|e^ -e-^| 2 

X £ 2T (Pn+l) | ^(Pl, " " ' ,Pn-l,Pn+Pn+l) ^ (4.6) 

«F, Vf~F)) = \j dMn+i(lWi) ^ 2T (Pn + i)0«(Pn + i)(e tri - e -^)( e « - e"*" + 1 ) 

X [F(pi +p n +l,P2 • • • ,Pn)F(p 1 ,p 2 +Pn+1, • " • ,Pn) + C.C.] (4.7) 

«F, ^ 2t F» = iy ^ n+1 ( Pn+1 ) ^ 2T ( Pn+ i)e K (p n+1 )(e^ -e-- 2 )(e-3 _ e -*r«) 

x [F(p l + p 2 , p 3 ■ ■ ■ , p n+ i)F(p 1 ,p 2 , p 3 + p 4 , ■ ■ ■ , Pn+l) + cc] (4.8) 
To check the combinatorics, we notice that the total number of terms are 
^ijn + l) ^ ^ + 2 ^ _ ^ + ^ {n- l)(n-2) \j = ^ n(n+ 1) \2 

the same as the total number of terms in (AF) 2 . The factors are obtained in the following 
way. Notice that in the formula of (AF) 2 we have to choose two indices. We first fix the 
special two indices in one F to be, say, (1, 2). This gives a factor n(n + l)/2. There is only 
one choice for the second index to be (1,2) and this gives the first factor for the diagonal 
term. There are 2(n — 1) choices to have either 1 or 2 and (n — l)(n — 2)/2 choices to have 
neither 1 nor 2. These give the last two factors. 

Notice that by the Schwarz inequality, the off-diagonal term is bounded by the diagonal 
term. For the purposes of upper bound we only have to estimate the diagonal term. Since 
the number of the off-diagonal terms are bigger than the diagonal terms by a factor of order 
n 2 , we have the upper bound 

J d/! n+ i(p n+ i) ^ 2r (p n+ i)0 K (p n+ i)|A + F(pi, • • • ,p n+ l)\ 2 

<Cn 4 ((F, K^F)) (4.9) 
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4.2 Preliminary Remarks 

Notice in the expression for we can integrate the variables p n — p n +i- So we make the 

change of variables and define some notations: 

u + =Pn+ Pn+1, U-=p n - p n +l, V^X = T n - r n+1 , \fcy = S n - S n+1 (4.10) 

Suppose at least one of \r n \, \r n +i\, \s n \, \s n+ i\ is not near or tt, say 

tt/100 < |r n | < 99tt/100 

Then we have 

) > oj{r n ) > C 

for some constant. Therefore, we can bound the kernel K (p ra +i) < and 

| e ir„ _ g-trn+i] = \ e i(r n +r n+ i) _ ^ < Cy/u{r n + r n+1 ) (4.11) 

After integrating p n — p n+ i, we change the variable u + = p n + p n +i to p n . Recall the 
normalization difference ((n + 1)!) _1 and (n!)™ 1 for dfi n+ i and dfi n . Thus we have 

J dfi n+1 (p n+1 ) {tt/100 < \r n \ < 99tt/100} 

e K (p n+1 )\e ir " - e ~ ir ^\ 2 | F(p u ■ ■ ■ ,p n -up n + p n+1 ) | 2 

<Cn- 1 J dfi n (p n ) u{r n ) | F( Pn ) | 2 (4.12) 

Since we are interested only in terms diverge as A — > oo, this term is negligible. Therefore, 
we shall assume that 

\r n \, \r n+1 \, \s n \, \s n+1 \ e [0, tt/100] U [99tt/100, tt] (4.13) 

We now divide the integration region according to \r n \, |r n+1 |, \s n \, \s n+ i\ in [0,7r/100] or 
[99tt/100, 7t]. There are sixteen disjoint regions and the final results are obtained by adding 
together the estimates from these sixteen disjoint regions. For simplicity, we shall consider 
only the region that all these variables are in the interval [0,7r/100]. The estimates in all 
other regions are the same. For example, suppose that r n+1 6 [997r/100, tt] and the other 
three variables belong to [0,7r/100]. Let p n+ \ = (tt, 0) + p n+ i and define 

G(p n ,p n+ i) = ^(Pn+l) 

Now we have |f n +i|> |§ n +i| ^ [0, 7r/100] and we can perform estimation on G instead on F. 
Notice that the factor e lTn — e~ lTn+1 may change by a factor of modular one. But we shall 
always take its absolute value in the estimates so that this factor disappears. 
Therefore, we now assume the following generality assumption 

Gi : \r n \, \r n+1 \, \s n \, \s n+1 \ E [0, tt/100] (4.14) 

This argument applies to all terms for the rest of this paper and we shall from now on consider 
only this case. The indices n, n + 1 are the two indices appear in F(pi, ■ ■ ■ ,p n -i,p n +p n +i)'i 
they may change depending on the variables we use in the future. Notice in this region, we 
have 

u(pj) ~ p) , j = n, n + 1, Lu(p n ± p n+1 ) ~ (p n ± pn+i) 2 (4.15) 

Since we concern only the order of magnitude, for the rest of the proof for Theorem \3.J\ in 
sections 4-6, we shall replace u(p) by p 2 whenever it is more convenient. 
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4.3 Upper Bound of the Diagonal Term: the Good Region 

The following Lemma is the main estimate on the diagonal term in the good region. 

Lemma 4.1 

((F, Kf 2T F)) < J dMPn) ^(r n )|log(A + u;(p n ))r- K/2 |F(p n )| 2 (4.16) 

Recall the change of variables ( |4.10p . We can bound the diagonal term from above as 



n-1 



.7 = 1 J 



(n+l)!^ E n-^. +u+=0j= 

tt/10 /-tt/10 

/ dxdy A + a 2 + b 2 + x 2 + y 2 + (a 2 + x 2 ) | log (A + a 2 + b 2 + x 2 + y 2 ) \ 

-tt/10 J -tt/10 L 



where 

b 2 = u>(s n _!) + u(e 2 ■ u + ), a 2 = w(r n _i) + io{e x ■ u + ). 
Clearly, we have 

g 2T ( Pn+1 ) C{x 2 + y 2 < C|logA|- 4 ^}{a 2 + b 2 < C|logA|- 4 -} 

We now replace u + by p n . Recall the normalization difference ((n + 1)!) _1 and (nl)^ 1 for 
dfi n+ i and d\i n . Thus we have the upper bound 

{{F, ^ 2T F))<-^— [ ^ n (p n ) co(r n )\F(p n )\ 2 {a 2 + b 2 <C\\og\\-^} 
(n + l)J 

x J J dxdy {x 2 + y 2 <C|logA|- 4r } 



x 



A + b 2 + y 2 + (a 2 + x 2 ){l + | log(A + a 2 + b 2 + x 2 + y 2 )\ K ) 



where 



b 2 = u>(s n ), a 2 = u>(r„) 
We need the following Lemma which will be used in several places later on. 

Lemma 4.2 Let r > 1 and 

Kl(a,b)= J J dxdy{x 2 + y 2 <\\ogX\- 2 ^} 

A + b 2 + y 2 + (a 2 + x 2 ){l + | log(A + a 2 + b 2 + x 2 + y 2 )\ K } 

Suppose that 



-i 



a 2 + b 2 < | log A | 



-2t 



(4.17) 
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Then for < k < 1 we have 



K T K (a,b)<C\ log(A + a 2 + 6 2 ) I 1 K ' 2 



On the other hand, if 



we have the lower bound 



a 2 + b 2 < |logA|- 4r , 



K T K (a,b) > C- l \ log(A + a 2 + 6 2 ) | W2 



We also have the trivial bound 




dxdy 



A + a 2 + b 2 + y 2 + x 2 



-i 



< C| log(A + a 2 + b 2 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



Proof: Clearly the trivial bound can be checked easily. We now prove the rest. Fix a 
constant m, 1 < m < r. Let 

G(x,y) = { (x,y) : \x\ <\\og\\ m \y\ < | log A| 2m |x| } 

and B be its complement. 

In the region B we drop (a 2 + x 2 )[l + \ log (A + a 2 + b 2 + y 2 )\ K } to have an upper bound. 
The angle integration of x, y gives a factor | log A| _m . Thus the contribution from this region 
is bounded by 

C|logA|-' n+1 < C 

In the region G we have 

log(A + a 2 + b 2 + y 2 )- 1 - log(l + | logA| 2m ) 
< log(A + a 2 + b 2 + x 2 + y 2 )" 1 
<log(A + a 2 + 6 2 + i/ 2 )- 1 



By assumptions ([4.17 ) or ( fi.19 ), a 2 + b 2 + x 2 + y 2 < 2| log A | 2r . Thus for r > m, we have 

log(A + a 2 + b 2 + y 2 )- 1 - log(l + | log A| 2m ) > Clog(A + a 2 + b 2 + y 2 y l 
for some constant depends on r, m. Therefore, we have 



C log(A + a 2 + b 2 + y 2 )- 1 < log(A + a 2 + b 2 + x 2 + y 2 )" 1 < log(A + a 2 + b 2 + y 



,2\-l 



(4.22) 



Upper Bound: We now replace log(A + a 2 + b 2 + x 2 + y 2 ) -1 by Clog(A + a 2 + 6 2 + y 2 ) -1 
and drop a 2 {l + | log(A + a 2 + b 2 + x 2 + y 2 ) | K } to have an upper bound for K^. Thus we can 
bound KJ.(a, b) by 



K T K {a,b) <C J J dxdy{x 2 + y 2 < |logA|' 2r } 

A + b 2 + y 2 + a 2 + x 2 \ log(A + a 2 + b 2 + y 2 )\ K ) 
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Change the variable by 



z = x 



\og(X + a 2 + b 2 + y 2 )\^ 2 



Hence 



z 2 < x 2 |log(A + a 2 + & 2 )| K 
Thus for x, y in the integration region we have 

z 2 + y 2 < |logA|- 2r |log(A + a 2 + 6 2 )| K < C 

We can bound K T K {a, b) by 




Kl(a,b)<C / / dzdy{z 2 + y 2 <C} 



A + a 2 + b 2 + y 2 + z 2 



log(\ + a 2 + b 2 + y 2 )r/ 2 



Denote by 



Since 



2 2,2 

p = z +y 



(4.23) 



log(A + a 2 + b 2 + y 2 )- 1 > log(A + a 2 + b 2 + p 2 )~\ 
we can bound the integration by 

C [ dp 2 (\ + a 2 + b 2 + p 2 Y l \ log(A + a 2 + b 2 + p 2 )|- K / 2 < C\ log(A 
Jo 



+ a 2 + 6 2X|l- re /2 



This proves the upper bound. 

Lower Bound: We now replace log(A + a 2 + b 2 + x 2 + y 2 )^ 1 by log(A + a 2 + b 2 + y 2 )^ 1 to 
have a lower bound for K T R . We change the variable to the same z and p as in the upper 
bound. We now restrict the angle 9(z,y) of the two dimensional vector (z,y) to be between 
7r/3 and 27r/3, i.e., 



tt/3 < 6(z,y) < 2tt/3 



(4.24) 



In this region, \z\ ~ \y\ ~ p. Denote by q = p 2 + a 2 + b 2 . We further restrict the integration 
to 



W = { 2(a 2 | log(A + a 2 + b 2 )\ K + b 2 ) < q < \ logA|- 2r /2 } 



(4.25) 



From the last restriction, we also have p 2 < \ logA| _2r /2. Since a 2 + b 2 satisfies ( |4. 19| ) and 
| a; | < \z\, the condition x 2 + y 2 < \ log A|~ 2t is satisfied. 
The integral is thus bounded below by 

/ dzdy {tt/3 < 9(z, y) < 2tt/3} W{q) \ log(A + a 2 + b 2 + p 2 ) \ ~ K ' 2 



A + b 2 + p 2 + a 2 \ log(A + a 2 + b 2 + p 2 ^ K 
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From the restriction on q, we have 

p 2 > a 2 \ log(A + a 2 + b 2 )\ K > a 2 \ log(A + a 2 + b 2 + p 2 )\ K 

Thus 



> (1/2)(A + a 2 + b 2 + p 2 ) 



2\-l 



A + b 2 + p 2 + a 2 \ log(A + a 2 + b 2 + p 2 )\ K 
The angle integration produces just some constant factor. Thus the integral is bigger than 

logA|- 2 V2 



C I ^(A + g)- 1 |log(A + g)r K / 2 

'2(a 2 | loga 2 | K +b 2 ) 



> c 



log(A + 2a 2 |loga 2 | K + 26 2 )| 1 K ' 2 - [2rlog | logA| + log2] 1 K ' 2 



Since a 2 + b 2 < \ log A| 4t we have 

log(A + 2a 2 |loga 2 r + 26 2 r 1 > (7r/2)| loglog A| 
Therefore, we have the bound 

|log(A + 2a 2 |loga 2 r + 26 2 )| W2 - [2rlog | logA| + log2] 1 ~ K/2 
> (1/20)| log(A + a 2 + b 2 )\ 1 - K/2 

We have thus proved the lower bound. 



(4.26) 



(4.27) 



Q.E.D. 



From this Lemma, we have proved Lemma |4.1| concerning the estimate in the good 
region. Observe that the main contribution of the p n — p n +i integration comes from the 
region \p n — p n+ i \ ^> \p n + p n +i \ + <-<^(Pn-i)- In fact, we have the following Lemma. 



Lemma 4.3 For any m > there is a constant C m such that 

<i/l n+ i(p n+ i) ||p„-p„+l| 2 < | log A| 2m [ \p n +Pn+l\ 2 + W(p n _l) ] j 



e ~ir„+i 12 



I F(P\, ■ ■ ■ ,Pn-l,Pn +Pn+l] 



A + w(p n+ i) 
< C m n~ x | loglog A| J dfi n (p n ) uj(r n ) \ F(p r , 



Proof: We have the bound 



d{p n -Pn+l) 



Ig^n _ g— WVi+1 12 



A + <^(p n +l) 

X {\Pn-Pn+l\ 2 < |logA| 2m [|p n +p n+1 | 2 + ^(p n _ 1 )]} 



< l0£ 



A + (l+|l0gA| 2m ) \Pn+Pn + l\ 2 + u(p 



n—1 



A + \p n + Pn+l | 2 + £*>(p n _l) 

Changing the variable p n + p n+ i — > p n , we have proved the Lemma 



< C m |loglogA| 



(4.28) 



Q.E.D. 



17 



Therefore, with a price of the term on the right side of 
general assumption (II) 



we can assume the following 



G 



11 



\Pn-Pn+l\ 2 > |logA| 2m [|p r , 



Pn+l\ 



+ w(p n _ 1 )] . 



(4.29) 



Under the assumptions ( |3.9|) (|3.10|) , the term on the right side of ( 4.28| ) is much smaller than 
the accuracy we need for Theorem O. Therefore this condition will be imposed for the rest 
of the paper. 



4.4 Upper Bound of the Diagonal Term: the Bad Region 

The contribution from the bad region can be decomposed into diagonal and off-diagonal 
terms. Again, we shall use the Schwarz inequality to bound the off-diagonal terms by the 
diagonal terms. Therefore, we have the bound 

\A+F( Pl ,--- , Pn+ l)\ 2 



J d/i n+1 (p n+ i) B 2r (p n+1 ) 



A + u>(p n+ i) 



/Igirn _ e ~ir n +i 12 
d/I n+1 (p n+1 ) B 2T (p n+1 )—— — | F(p u ■ ■ ■ ,p n -l,Pn+Pn+l) f 
A + U)(J) n+1 ) 

Again the variable p n — p n +i does not appear in F and we can perform the integration. 
We subdivide £> 2r (p n+ i) into 

B 2T (p n+1 )Bn T (Pn-l,Pn+Pn+l) U B 2t (p n +l) (Pn-1 , Pn + Pn+1 ) ■ 

In the first case, we drop the characteristic function B 2r (p n+ i) to have an upper bound. We 
now use the trivial bound (|4.21|) to estimate the integration of the variable p n — p n +i by 

2n 4 J dfJL n+1 (p n+1 ) B 4T (p n - 1 ,p n + p n+1 ) 

| e ir„ _ e -ir„ +1 |2 
X — r- 7 r~ F(pi, ■ ■ ■ ,Pn-l,Pn+Pn+l) 

X + w(p n+ i) 

< Cn 3 j d^ n {p n ) B iT {p n )uo{r n )\ log(A + w(p n ))| | F( Pl , ■ ■ ■ ,p n ) \ 2 

< Cn 3 \ log log A | / dyL n {p n ) u(r n ) \ F(p u ■ ■ ■ ,p n ) \ 2 (4.30) 



Here we have used the change of the normalization between dfi n+ i and d/j, n . 

We now estimate the region i3 2T (p n+ i)^ 4T (p ra _i, p n + p n+ i) which is the transition from 
the bad set to good set. In this region, 

bn-Pn + l| 2 >C|logAr 4T 

The contribution is bounded by 

2n 4 y d/i n+1 (p n+ i) {\p n -p n+1 \ 2 > C|logA|" 4T | 

X — r- 7 r~ I F(pi, ■ ■ ■ ,Pn-l,Pn+Pn+l) f 

X + u>(p n+1 ) 

<Cn 3 |loglogA| J dfi n {p n ) u{r n ) \ F{p n ) \ 2 (4.31) 
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Combining the estimates ( [4.30 ) and ( [4.31 ), we can bound the contribution from the bad 
region by 



d/in+l(Pn+l) B T (Pn+l) 



\A + F(p U --- ,p n +l)\ 

A + u>(p n+ i) 



< Cn i \ log log A | / dp, n (p n ) uj(r n ) \ F(p r , 



Together with the estimate on the good region, Lemma O] , we can bound the right side of 
(0) by 



^n+l(Pjt+ly 



\A + F( Pl ,--- ,p n+1 )\ : 



A + a;(p n+1 ) + V+ 2 " r +1 (p n+1 ) 
<Cn 3 J dfinipn) u(r n )\ log(A + u,(p n ))| 1 - K / 2 |F(p„)| 2 

+ Cn 3 \ log log A| J d/i„(p„) uj(r n ) | F(p n ) | 2 

Under the condition ( |3.10| ), it is easy to check that for symmetric function F the right side 
of the last equation is bounded above by | loglog A| 2 U~ T (p n ). This proves the upper bound 
for Theorem I3TT . 



5 Lower Bound: The Diagonal terms 

By definition, we have 



({FA^X-S^ + ^V^r'AF}} 



\A+F(p U ■■■ ,p n 



+1, 



A + w(p n+1 )+ 7 - 1 U^ 1 (Pn+i; 

Since 7 < 1 and A + uj > 0, the integral is bigger than 

\A+F{ Pl ,--- , Pn +l)\ 2 



(5.1) 



7 / dfi n+1 (p n+1 ] 



A + U>( Pn+1 ) +U«+ 1 (Pn+l) 



Divide the integral into p n+ i G Q T and p n +i G £> T . In the bad set B T , we bound the 
integral in this region from below by zero. In the good set, we have 



U^(p n+1 ) = u>(r n+1 )| log(A + o;(p n+1 ))r, p n+1 e Q 1 



Thus 



((^(A-^ + U^ 1 )- 1 ^)) 



> / d/i n+ i(p n+ i) G T (p n +i)- 



\A + F{ Pl ,--- , Pn+1 )\ 2 
'A + u;(p„, +1 ) + U^+ 1 (Pn + i) 

> / d/j n+1 (p n+1 ) Q T (p n+1 )Q K (p n+1 )\A + F(p 1 , ■ ■ ■ ,p„+i)| 2 



(5.2) 
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where K (p n +i) is defined in ( |4.4j ). We now decompose the last term into diagonal and 
off-diagonal terms: 

+ n ( n -l)( n + l)«ir, f^F» + " ( "" 1)( " 4 " 2)( " +1) ((F, Vfn (5.3) 
where these operators are defined in ( [4.6| )-( fX~8| ). 

5.1 Lower Bound on the Diagonal Terms 

The main estimate on the lower bound of the diagonal term ( |4.6[ ) is the following Lemma. 
Define 

F<T(p n ) = F(Pn)G 2T (Pn), F^( Pn ) = F(p„)i^(p„) (5.4) 



Lemma 5.1 Recall k,t and n satisfy the assumptions ( |3.9| ) and ( |3.f0| ). Then the diagonal 
term is bounded below by 

({F, K^F)) > Cn- l ((F^, u{r n )\\og{\ + u{ Vn ))\ l -^F^)) (5.5) 



Proof: Recall the assumptions ( |4.14j ), ( |4.29|) and the change the variables 



u + = Pn+ Pn+1, U- =p n - p n+ i, V2x = T n - T n+1 , \^2y = S n - S n+1 

b 2 = u(e 2 ■ u + ) + u;(s n _i), a 2 = Lo{e x ■ u + ) + w(r n _i) (5.6) 



Thus we can bound the diagonal term from below by 

v ii 

n-l 



[F, F\ 



> , ~ rvj / I [dPj I du + uj(e 1 -u + )g T (p n+1 )\F(p 1 ,--- ,p n _ u u + )\ 

X + a 2 + b 2 + x 2 + y 2 + (a 2 + x 2 )\ log(A + a 2 + b 2 + x 2 + y 2 )^' 



2 



x J J dxdy 

We now impose the condition x 2 + y 2 < | log A| _2r /2 to have a lower bound. Since 
Q 2t (Pi, ■ ■ ■ ,Pn-i, u+) { x 2 + y 2 < I log \\- 2t /2 } C g T (p n+ i), 

we can replace G T (p n +i) by { x 2 + y 2 < | log A| _2t /2 } and F by Fg T to have a lower bound. 
The lemma now follows from the lower bound of Lemma [0. Q.E.D. 
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6 Off-diagonal terms 



Our goal in this section is to prove the following estimate on the off-diagonal terms. 



Lemma 6.1 Recall that k,t and n satisfy the assumptions ( |3.9|) and (|3.10|) . The first and 
second off-diagonal terms are bounded by 



KCvr^loglogX^ 1 ' 2 J df, n {p n ) u(r n ) |F^( Pn )| 2 

+ Cn' 5 J dfi n (p n ) ^(r n )|log(A + a;(p n )| 1 -' 1 / 2 |F^(p n ) 



Proof: The first off-diagonal term is bounded by 



<C j dfM n+1 (p n+1 ) Q T (p n+1 )Q K (p n+1 )\(e^ - + 1 )(e^ - e"™)| 

X F(p 1 +p n+ i,p 2 ■ ■ ■ ,Pn)F{p 1 ,p 2 +Pn+l, ■■■ ,Pn) 

By definition F = Fg T + Fj r . Thus the last term is equal to 



C j dfi n+1 (p n+1 ) £ r (Pn + i)e K ( Pn+1 )|(e^ - e ~ tr ^)(e" 

X {Fg T + F% T ) (p! + Pn+1 , p 2 , ■■■ ,p n ) 
X (F$ + F£) (p 2 + Pn+1 , Pl ,p 3 , ■■■ ,p n ) 

From the Schwarz inequality, the cross term is bounded by 



«f 2 _ -ir n+1 



C J df, n+1 (p n+1 ) £ T ( Pn+1 )e K ( Pn+1 )|(e iri - e - ir ^)(e ir * - 

x Fg T (p 1 +p n+1 ,p 2 ,--- ,p n )Fg T (p 2 +p n+1 ,p 1} p 3} --- ,p n ) 
<C5 J ^„ + i(p„ + i)^ r (p„ + i)e ft (p n+1 )|(e iri -e- ir "+ 1 )| 2 

2 

X Fg T (p 1 +p n+1 ,p 2 ,--- ,p n ) 

+ CS- 1 J dfx n+1 (p n+1 ) g T (p n+1 )e K (p n+1 )\(e ir2 - e - ir » +1 )| 2 

2 

x Fg T (p 2 +p n+1 ,p 1 ,p 3 ,--- ,p n ) 
We first bound the last term. Clearly, in the region 

Q T (Pn+l) S 2T {p 2 +Pn+l,Pl,P3, ■■■ ,Pn} 



(6-1) 



(6.2) 
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we have 

\P2 ~Pn+l\ 2 < |logA| 4r [|p 2 +p n +l\ 2 +W(pi) + Uj(p 3 ) H w(p„) ] 

Thus we can apply Lemma L3. Let 5 = | log log A| _1//2 . We can bound the last term in ( 6.2 
by 



Cn- 1 ! loglogAr +1 / 2 /" d// n (p n ) w(r n ) |i^(p, 



)| 



(6.3) 



The first term on the right side of 
Using Lemma |4.1| , we bound it by 



can 



be bounded as in the section of upper bound. 



Cn- l \ loglog A)- 1 / 2 J d» n (p n ) u(r n )\ log(A + p^l W2 ^(Pr. 



(6.4) 



The contribution from the term with Fg T Fg T can be estimated similarly. Finally, we 
consider the contribution from Fg T Fg T . To estimate this term, we need the following Lemma 
which will be proved in the next section. 



Lemma 6.2 Recall that k,t and n satisfy the assumptions ( |3.9|) and ( p.lOj ). Then we have 
the following two estimates: 

Qi = J <Wi(Pn+x) Q K (Pn+i)\(e iri - e-^)(e^ - e^)\ 

Fg T (pi+P3,P2,P4,--- ,Pn+l)Fg T (p 2 +P3,PuP4,--- ,Pn+l) 

< Cn~ 5 J d/i n (p n ) cu(r n )|log(A + a;(p n ))| 1 - K / 2 |F 2 ^(p„)| 2 (6.5) 



Q 2 =J d/Wi(Pn+i) e K (p„ +1 )|(e iri -e- ir2 )(e ir3 -e— 4 )| 

Fg T (Pl+P2,P3,P4,--- ,Pn+l)Fg T (p 3 +P4,Pl,P2,P5,--- ,Pn+l) 

< Cn' 5 J dfx n {p n ) w(r n )|log(A + w( Pn ))| 1 - K / 2 |F 2 ^(p„)| 2 



(6.6) 



We now collect all our efforts. The cross terms are bounded by ( |6.3|) and ( |6.4|) . The 
contribution from Fg T Fg T can be estimated similarly. Finally the contribution from Fg T Fg T 
is bounded by the last Lemma. Thus we have proved the estimate on in Lemma |6.1| . 

The estimate on can be proved in a similar way by using instead the equation (|6.6|) . 

This proves Lemma |6.1| . 
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6.1 Proof of the Lower Bound 

Recall the condition ( |3.10|) on the size of n. Combining the lower bound on the diagonal 



term in Lemma |5J] and the estimate on the off-diagonal terms in Lemma |6.1| , we have 



n 3 |«F, ^F\ 



n 4 \((F, V^F] 



>Cn 2 J dfi n ( Pn ) u;(r n )|log(A + w( Pn ))| 1 - K / 2 |F^( Pn )| 2 
- Cn 4 |log log A| 1+1 / 2 y dfj, n (p n ) u{r n ) \F^( Pn )\ 2 

_ dn' 1 + n 3 \ log log A)" 1 / 2 ] / d// n (p n )u;(r n )| log(A + ^(p^))! 1 — / 2 |F^(p n )f 



The last term can be absorbed into the first term on the right side with a change of constant. 
The middle term on the right side gives the estimate on the bad set. This proves the lower 



bound for Theorem 3.1 



6.2 Proof of Lemma |6.2| 

We first bound Q\. Consider the two cases. 

Case 1. Some pi, i = 1, 2, 3 dominates, say, we have 

|pi|>2(|p 2 | + |p 3 |) 

Then |pi — p 3 1 < 4|pi +p 3 |. From the Schwarz inequality 

J c^„+i(p„ + i) {|pi -p 3 | < 4|pi + p 3 |} 



x Q K (p n+1 ) (e^-e- ir3 )(e^-e 



X Fq T (V\ +P3,P2,P4, ••• ,Pn+l)Fg T {p 2 +P3,Pl,P4,--- , Pn+l] 
< 5~ l J dfji n+1 (p n+1 ) {|pi -p 3 | < 4|p! +p 3 |} 

x e ft (p„ +1 )|e^-e-^| 2 |F 2 -(p 1 +p 3 ,P2,P4,' 
+ 5 j dfx n+l (p n+1 ) e K (p„ +1 )|e ir2 -e-^| 2 

x \ F g T (P2 +P3,Pi,Pa, ■ ■ ■ ,P«+i)| 2 



(6.7) 



The last term on the right side of ( |6.7|) can be bounded using Lemma |4.1] . To estimate 
the first term, we drop u>(r n+1 )\ log(A + o;(p ri+1 ))| K in K (p n+ i) and integrate pi — p 3 . The 
integration can be estimated easily by 



J d(pi-p 3 ) {IPi-PsI <4|pi+p 3 |} |©«(Pn+i)| <C . 



(6.8) 



We now choose 5 = n 5 and use 



n 



10 



<|loglogA| 1 / 2 <C|log(A + u;(p n )) 



l-ft/2 
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if u>(p n ) < | log A | 4r and < k < 1. The left side of ( |6.7[ ) is thus bounded above by 

Cn- 5 J d/x n (p n )w(r n )| log(A + u;(p n ))| W2 |F 2r (p n )| 2 

Here we have changed variables so that the variable of the function Fg T is of the standard 
form. 

Case 2: \p x \ ~ \p 2 \ ~ |p 3 |- 

In this case, we have |pi — p 3 1 < 16 1 j?2 1 • Similar arguments prove the same bound in this 



region. This proves 

We now estimate Q 2 - We can assume without loss of generality that 

<*>(pi - Pi) < u{p z -pi) 
Again, we bound it by the Schwarz inequality to have 

Q 2 <6- 1 J dfi n+1 ( Pn+1 ) e K (p n+1 )\e^-e'^\ 2 

X {^(Pi ~P 2 ) < ^(P3 ~ P4)}\Fg T (pi + p 2 ;P3,P4, " " ' ,Pn+l)\ 

+ 5 J dfi n+1 (p n+1 ) & K (p n+1 )\e ir3 -e~ ir *\ 2 

X \Fg T {P3+P4,Pl,P2,P5,--- ,Pn+l)\ 2 

Both terms can be estimated by similar arguments used for Q\. So we obtain ( |6.6| ). Q.E.D. 



|2 



7 Conclusions 

From the main estimate Theorem we need the relation 

K n _\ = 1 — K n /2. 

To satisfy this relation, for any large integers N fixed, we let 

K n = 2/3 + (-l) n 2- 2JV+n /3, n = l,--- ,2N + 1. (7.1) 
A few terms are given explicitly in the following: 

^2iv+i = 2/3 - 2/3 = 0, k 2N = 2/3 + 1/3 = 1, « 2Ar -i = 2/3 - 1/6, 

/t 2iV _2 = 2/3 + 1/12, • ■ ■ , «2 = 2/3 + 2- 2Af+2 /3 
We first apply Theorem |3J] to have 

A* D^ +1 A + <C\loglog\f\Jl» 
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In order to satisfy the condition 7 < | log log A |~ 3 later on, we now replace | log log A | 2 on the 
right side by | log log A| 3 to have a further upper bound. Now we apply the lower bound part 
of Theorem |3.1| to have 

A%\d2n + A* + D^ +1 A + } A + > C\ loglogA^V 2 ^ 

We can repeat this procedure until we have 

A* + (D 3 + ■ • ■y 1 A + < C\ loglog Ar+ 4 l4 iT 

Thus we have 

((w, [Dz + AKDz + .-.y'A+Y'w)) > ((w,[D 2 + C\ loglog Al 2 ^ 4 ^^]"^)) 
The Fourier transform of w is 

w(pi,p 2 ) = e~ lT2 
Since p\ + p 2 = under the measure dfi 2 , we have 

{{w,[D 2 + C\ loglog A| 2JV+4 U 2 2 J-V ; )) 



=1 j dpi {A + 2^(p 1 ) + C|loglogA| 2Ar+4 U 2 2ir (p 1 ,-j9 1 )]|^ 1 



The last integration is the same as the right side of ( |5.1|) with n = 1 and A + F replaced 
by one. Following similar argument, we have 



^ T 2 (o,o) 



((w, [D 2 + C\ loglog Ar+Xj-M) > C\ loglog \\- 2N - 1 * 
where K^ 2 (0,0) is defined in Lemma L2. From ( 4.20| ), we have 

^(0,0) > |logA| Kl , k x = 2/3 - 2- 2N+1 /3 

Thus we have 

((w, [D 2 + C\ loglog AI^U 2 ^]- 1 ^)) > I loglog A|- 27V - 4 | logAr 
Therefore, we have the lower bound 



D 2 + A* + (D 3 + ---) X A + 
I log log A 



((w 

> I log A| 2/3 exp 



By choosing 



2 2N-1 3 

N = a I log log log A 



w)) 

- (2N + 4)\ log log log A I 



with a large enough, together with Lemma |2.1| we have proved the lower bound. 



Instead of ( |iM|) , we can choose 

n n = 2/3 - (-l) n 2- 2N+n+1 /3, n = l, 



,2N. 
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Explicit examples are 

k 2N = 2/3 - 2/3 = 0, ft 2A r-i = 2/3 + 1/3 = 1, k 2N _ 2 = 2/3 - 1/6, 

k 2 jv-3 = 2/3 + 1/12, • • ■ , k 2 = 2/3 - 2- 2Af+3 /3 
With this choice of K n , similar argument proves the upper bound. This concludes Theorem 

m 
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